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By Robert E. Davidson 

SUMMARY 


The Bernoulli equation for nonsteady flow was used to derive 
formulas for propeller lift and thrust distribution in terms of wake- 
survey measurements of stagnation pressure rise through the propeller. 
Lift distributions and over-all values of thrust from wake survey were 
compared with direct measurements of these quantities at stream Mach 
numbers up to O. 65 . 


INTRODUCTION 


In reference 1, Lock and Yeatman derive the formula for propeller 
thrust in terms of wake-survey measurements. Their formula is the same 
as is obtained from momentum considerations, but the method of derivation 
is believed to give greater insight into the mechanics of the wake survey. 
The analysis of Lock and Yeatman, although complicated at the start by 
the introduction of the velocity potential, develops the thrust formula 
in a logical manner for a finite number of blades and periodic flow. 

From this development it appears that thrust distribution is not the 
only thing which the wake survey is capable of measuring. For instance, 
it can be deduced that lift distribution is given more naturally by the 
wake survey than thrust. Reference 1 shows that the derivation of the 
thrust formula actually involves neglecting the induced angle of attack; 
in fact, once this assumption is made, the wake survey can just as 
reasonably be used to measure torque as thrust. 

In this paper it is shown that the wake -survey formulas may be 
modified for compressible flow without very complicating changes in 
form. As in reference 1, the total -pres sure survey tube is recognized 
as measuring the time rate of change of the velocity potential. If 
advantage is taken of the fact that pressure changes are small in the 
flow about a propeller, the wake -survey formulas for compressible flow 
remain simple in form. 
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SYMBOLS 


number of blades 


blade chord 


power coefficient 



thrust coefficient 



section lift coefficient 
section normal -force coefficient 
section power coefficient (dCp/dx) 
section thrust coefficient (dCT/dx) 
propeller diameter 



(See equation (9)*) 


— (See equation (9)*) 


advance ratio (v/nD) 


advance Mach number 


helical section Mach number 



rotational speed, revolutions in unit time 
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p static pressure 

p^ stagnation pressure 

r radius from the rotational axis to a point in the flow 

T thrust 

t time 

V advance velocity 

W resultant velocity of a propeller section 

Wq velocity vector (see fig. l) 

w fluid velocity 


X 


radius divided hy propeller radius 



r 

P 

pt 

0 

O) 

T 

7 

Apt 

Subscripts : 


circulation 
mass density 
stagnation density 
velocity potential 

rotational speedy radians in unit time 

time between blades passing a survey tube (2Jt/Bo)) 

ratio of specific heats 

stagnation pressure rise through propeller 


1 denotes free stream condition far ahead of propeller or at 

some place where flow is practically steady 

I ) curly brackets j | are used to show that a time average is 

’ ^ taken ^ ^ 
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DERIVATIONS 


The flow field about a propeller, as for a wing, is irrotational 
with the exception of the thin surfaces of rotational flow which trail 
from the blades as a result of viscosity. When the effect of the blade 
wakes on a survey tube is neglected by assuming the flow to be all 
irrotational, the velocity, the velocity potential, and the state 
functions of the gas must satisfy: 

M=E+iw^+E - F(t) (1) 

St P 2 

where 

E = - Jv d(l/p) (2) 

and the coordinate system is fixed to the survey rake. 

Observe that for the propeller wake survey, unlike the wing, the 
irrotational flow outside the vortex sheets is surveyed and the measur- 
able forces are those arising primarily from circulation rather than 
viscosity. 

Equation (l) is applicable to the compressible flow past a propeller 
with no restriction on the magnitude of the velocities except that the 
velocity potential 0 must exist. This equation appears on page 20 
of reference 2. 


For a compressible flow the following relation still holds: 

r = I c^bw (3) 

where F is the circulation at a particular radius and is equal to the 
jump in 0 across a trailing vortex sheet at that radius in the wake. 

The relation between pressure and density is assumed to be as 
follows: 


— = Constant 
7 


( 4 ) 


which is known to be a good assumption for practical problems in the 
aerodynamics of compressible flows. 
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In compressible flow the change in velocity potential between the 
departure and arrival of a vortex sheet at a total -pres sure tube of the 
survey rake (see fig. l) is still equal to the jump across a vortex 
sheet which is equal to the circulation, all at a particular radius. 
Therefore, the reasoning of reference 1 with equations (l), (3)^ and 
(k) gives for conditions at a total-pressure tube behind the propeller 



whereas, far upstream. 


and F(t) = 0 because it is zero far ahead where the flow is steady. 
Equation (4) substituted in equation (2) gives 


where Pq, Pq represent some initial state. Subtracting equation (6) 
from equation (5) and using the expression for E gives 


The stagnation density behind the propeller is replaced by using 

equation to give 


When this equation is substituted in equation (7), the latter takes the 
form: 



(6) 




( 7 ) 




(8) 
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Although equation (8) gives the circulation in terms of total- 
pressure measurements alone^ it is desirable to be able to use the 
reading of an ordinary liquid manometer which is usually assumed to 

give /PtK rather than the averaged quantity in equation (8). A binomial 
7-1 

expansion of pt ^ will avoid the difficulty. Only one term of the 
series need be retained in order to be in keeping with the error in the 
experimental measurement of p^. Retaining this term leads to a formula 
which permits easy computation of the lift from total-pressure measure- 
ments with a liquid manometer. 

Let: 


pt = Pti + ^Pt 


h 


APt 

Ptl 




( 9 ) 


m = 



u 


then equation (8) becomes 


r 

T 


1 Hi 




and with the expansion of 


(l + h)^ this becomes the alternating series 


£ = i Eii 

T m 

In trying to see the relative magnitude of the terms, it is helpful 
to make the exponent in the second term of equation (lO) appear on the 
outside rather than the inside of the curly brackets. For this purpose, 
some assumption is necessary as to the form of h. 

Measurements with instantaneous total -pres sure measuring devices 
indicate that the total -pres sure variation with time between vortex . 
sheets is approximately parabolic. A typical measurement of this kind 
is shown in figure 2(a). The pertinent parts of the record are shown 
in figure 2(b). As would be expected, the total -pressure increment Ap^ 
apparently never comes down to zero at the midpoints. However, a 


m 




21 \ i 


( 10 ) 
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conservative simplifying assumption is that this does happen as shown 
in figure 2(c), that is, that the origin of the parabola is at Ap-j. = 0. 

With the assumption that h varies in time as represented in figure 2(c), 
the quantities (h|2 and | h2| are no longer independent, in fact 

!h2| = 

^ ^ 80 i * 

as some simple integrations will show. Equation (lO) then becomes 


p 1 PtlF < i m(m-l) l^^hip ~I /-.-.s 


For the values 


m = = 0.286 


p^2. ~ 2000 pounds per square foot 
jAp^l = 20 pounds per square foot 
|h| = 0.01 

which give a conservatively heavy disk loading, equation (ll) becomes 


r 

T 


[0.00286 - 0.000018 + 

m 



Therefore, the error when retaining only the first term in equation (ll) 
is less than 


0.000018 

0.00286 


= 0.0063 or 0.6 percent 


It follows that, to a good approximation. 
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or 


(i; 

T Ptl 

Equation (12) gives the circulation distribution from easily 
measured quantities. The bracketed term is obtained directly from a 
liquid manometer, whereas is the stagnation density in the 

undisturbed stream. 


Lift coefficient.- The lift distribution is easily obtained from 

equation (12). The time for one blade to succeed the preceding one is 
the same as the time between passages of the vortex sheets at a total- 
pressure tube} therefore, 


T 



(13) 


Then equations (3) and (13) substituted in equation (12) give 


ic^ bWBoD 
2n 


Ptl 


with which the lift coefficient may he calculated from stagnation 
measurements . 


With the excellent approximation that W = Wq (see fig. 1), the 
last equation may he solved for c^ 


Cl = 


b(>,/d) 


_i^ 

2 


(14) 


which is the formula with which c^ was calculated in figure 3* This 
formula will give slightly low values of c^ depending upon the extent 
to which profile drag influences the readings of the total -pres sure 
tubes. Because the radical expresses the section velocity at the blade, 
X must be at the blade with no correction factor concerning the slip- 
stream contraction. 

Thrust.- By projection of the lift in the forward direction, the 
thrust may be found from equation (l4). The induced angle (see 


2Z 
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fig. 1) and the drag are neglected. The formula for the local thrust 
coefficient in compressible flow is then 


H = 






(15) 


which is the same as the formula at the bottom of page 6 in reference 3. 

Integration of equation (15) gives the over-all thrust coefficient 
Cip shown in figure h. 

For incompressible flow, equation (15) reduces to the well-known 
formula for thrust 


— . «xD|APt| ( 16 ) 

which is obtained from the momentum theorem for incompressible flow. 

The substitutions 

dT = p3_n2DVc,p 


and 


dr = ^ D dx 


give equation (l6) in notation comparable with equation ( 15 ): 


_ i j2 


1 




(17) 


Equation (17) is for incompressible flow. Equations (15) and (17) are 
the same for incompressible flow because, in this case, the densities 
and are equal and all other terms are the same. Thus, the 

only change needed to make the established thrust formula ( 17 ) apply to 
compressible flow is the substitution of stagnation density for 

stream density p^. 


Formula (l6) is known to neglect the rotation of the slipstream 
when derived from momentum thoery. Formula (l6) neglects the slipstream 
rotation and formula ( 15 ) neglects the induced angle of attack. 
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COMPARISON OF WAKE -SURVEY RESULTS WITH PRESSURE 
AND OVER-ALL FORCE MEASUREMENTS 


An experimental check of the wake -survey formula (l4) for lift 
coefficient is possible. Propeller tests in which normal-force coeffi- 
cients were obtained at various stations on a 10 -foot -diameter propeller 
have been made in the Langley l6-foot high-speed tunnel. The normal - 
force coefficients are computed from direct measurements of static- 
pressure distribution on sections of the rotating propeller. Wake 
surveys were also made during the same tests. 

Lift coefficients from wake -survey measurements and formula (l4) 
are compared with normal-force coefficients from direct measurements 
of static pressures in figure 3- The tests from which the data were 
taken are reported in reference 4. In figure 3^ solid lines are faired 
through points plotted from wake -survey data. Plotted points with 
crosses through them denote normal-force coefficients obtained from 
direct static-pressure measurements on the blades. The points with 
crosses should fall on the faired lines except for the small difference 
between c^ and c^. The dashed lines are faired in through radial 
stations where a given sectional Mach number prevails. Disagreement 
is seen to come in around Mx = 1.0. 

In figure 4, thrust from wake-survey formula (15) is compared 
with direct thrust measurements. The comparison is not as thorough as 
that for the lift coefficient because only the integrated value of the 
wake -survey thrust is shown. The tests from which the data were taken 
are reported in reference 5* The agreement is fair up to the highest 
advance Mach number tested. 


CONCLUDING REMARKS 


Wake-survey formulas were presented which are applicable to 
compressible flow and convenient for computation if a slight degree of 
approximation is accepted. The usefulness of the wake survey has always 
been in the ease with which force distribution may be obtained from 
easily measured quantities. Therefore, a small sacrifice in accuracy 
was considered worthwhile in view of the resulting simplification of 
computations . 

As long as a velocity potential exists no Mach number restrictions 
are required. However, if for some reason the flow becomes seriously 
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rotational^ the formulas presented are in error by some indefinite 
amount . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Fields Va. 
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Figure 1.- Schematic diagram of propeller wake and survey setup, 
diagram of velocities at the blade. 
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(a) 



(b) 



APt 

Figure 2.- Actual and assumed time history of h = . 

Ptl 


Normal- force coefficient c Normal-force coefficient 
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(a) llto rpm; Pq .731 ~ 



(b) 1350 rpnij pQ , 75R ” ^5*^’ 


Figure 3.- Comparison of lift coefficient c^ from wake svirvey with 
pressure distribution normal-force coefficient c^. 



Normal- force coefficient, Oj^ Normal force coefficient 
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(c) 1500 rpm; = 45°. 


fe 


O 






(d) l600 rpm; Pq^.^^ = 45°. 


FigTjre 3.- Continued. 
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Figure 3-- Concluded 


ThruBt coefficient 


CO 

tSl 




Advance ratlo^ J 


(a) llkO rpm; Po.75R = 30°. 


(b) 2160 rpm; Po.75R 30°. 


Figure 4.- Comparison of thrust coefficient from wake survey with that 

from over-all force test. 
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Thriuat coefficient 




(c) llto rpm; = ^5°. 


(<i) Mq 


0.56; 


^0.75R 


45°. 


Figure 4.- Continued. 
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Advance ratlo^ J 


(e) Mq = 0 . 60 ; = 45°. 


Figure 4 


Thrust coefficient 
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(f) Mq = 0.65; = 45°. 
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Concluded 
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